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Abstract 
We compute the Casimir energy of the photon field in a static de Sitter space and find it to be 
proportional to the size of the horizon, the same form of the holographic dark energy. We suggest 


to make metamaterials to mimic static de Sitter space in laboratory and measure the predicted 


Casimir energy. 
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I. INTRODUCTION 


The dark energy problem is one of the central unsolved problems in fundamental physics 
and cosmology. Many phenomenological models have been proposed since the discovery 
of the accelerated expanding universe hh, including the so-called holographic dark energy 
model [2]. 

The value of the observed dark energy is comparable to the critical energy density which 
in the Planck unit is proportional to the Hubble constant squared. Motivated by this fact, 
the holographic dark energy model presumes that the energy density is inversely proportional 
to the square of the size of some cosmological horizon, and the choice of the horizon, for the 
model to work, is proposed in |2] to be the event horizon. In a pure de Sitter space in which 
dark energy is a constant, one can choose a static coordinate system, thus the total energy 
of dark energy can be defined. The total energy in the Planck unit is proportional to the 
size of the horizon, unlike the usual result for the Casimir energy in the finite cavity in a flat 
spacetime, the latter is inversely proportional to the size of the cavity. Nevertheless, there 
has been no attempt to computing the Casimir energy in a static de Sitter space seriously, 
to the best of our knowledge. 

In this paper, we take the assumption that the UV divergent zero-point energy is somehow 
regularized to vanish due to some unknown mechanism, therefore the observed dark energy 
is solely due to the residual zero-point energy which is known to be the Casimir energy. 
However, we are faced with the challenge that why the Casimir energy is not inversely 
proportional to the size of the cavity. 

Motivated by the recent development in electromagnetic metamaterials, we will first try 
to convert the de Sitter space to a cavity of optical metamaterial, this is doable at least for 
the photon field. We shall find that the cavity is drastically different from a usual one in 
that the permittivity and permeability parameters both have divergent components tangent 
to the boundary, indeed these components are divergent near the boundary which is derived 
from the horizon of the de Sitter space. This fact encourages us to guess that the Casimir 
energy as a function of the size of the cavity behaves differently from that in a normal cavity. 
We carry out the calculation for the photon field and find that indeed after regularization 
there is still a divergent piece proportional to the size of the cavity, with the divergence 


cut-off by a UV scale. This UV scale in the metamaterials can be some intrinsic microscopic 


scale, and becomes the Planck scale when transformed back to the de Sitter space. Thus, 
we find a piece of theoretical evidence for the holographic dark energy model. We propose 
to do some laboratory experiment to verify our result about the Casimir energy. This kind 
of experiment is interesting in two aspects. First, it will mimic cosmology, second, it will 
detect a unusual Casimir energy inversely proportional to some microscopic cut-off. 

There are two different designs of metamaterials, the first is based on the static de Sitter 
space with the usual coordinate system, the second is based on the static de Sitter where 
the radial coordinate is replaced by the proper distance. The Casimir energy in the second 
metamaterial assumes the precisely same form as in the de Sitter space, which is proportional 
to the size of the material thus is enhanced to be easily detectable. 

In the next section, we present the Maxwell equations in a curved spacetime. We discuss 
how to transform the de Sitter space into a cavity of metamaterial in sect.3 and compute 
the Casimir energy in sect.4, and discuss its form in a metamaterial filled cavity in sect.5. 


We conclude in sect.6. 


Il. MAXWELL EQUATION IN STATIC CURVED SPACE TIME 


We start with Maxwell equations in a curved spacetime 
pF =0, 0 H” = 0, (1) 


where the first equation is the Bianchi identity satisfied by field strength F», and H™ is 
defined by 
H" = /—gF” = y=gg"° g”? Fag. (2) 


The electric field Æ; and magnetic field H; are related to Fy as 
B,= Fy, Hy = -2E (3) 


In above equations, the Levi-Civita symbol ¢* is +1 for all even permutations of €!?°. It 
should be noted that E; and H; are spatial vectors. 
Especially, when the space time is static, we can express Eq. in a spatial covariant 


form with respect to optical metric y? defined by 


a” = —goog"’, (4) 


where Ņ* is the usually called optic metric. 


At this stage, Maxwell equations Eq. (I) can be rewritten as 


Vie =), V:H’ =0, (5) 
etik l etik 
oE’ — vate = 0, oH’ + mo = 0, (6) 


where V; denotes the covariant derivative with respect to y’’, and all the indices are raised 
and lowered by y. 


Eq. (6) can be rewritten in the following compact form 


E E 
iO; =f 7 (7) 
H H 
where the operator L is given by 
0 y ðk 
L=i ml va (8) 
Vij Ok 0 
If we define 
E 
o= (9) 
H 


the operator L is Hermitian in terms of the following inner product 


< Qı, Q2 >= | Ea aoto = | PaE Bs + Huu” Haj), (10) 
where 
S = fay", y = det (ij). (11) 


We will see that € and u have a physical interpretation in terms of metamaterials in the next 
section. 

To make L Hermitian under the inner product (10) (which is the energy functional as it 
should be for ¢; = $2), we need to specify boundary conditions if there is a boundary. We 


find the following boundary conditions 
Viale x H» = Ay x E>)'nilam = 0 (12) 


where n; is the normal vector of boundary and (E x H\ = Ae * FH, Then the hermiticity 
of L is guaranteed 


< Qı, Log >=< Loi, Q2 >. (13) 
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By acting ô, and L again on the two sides of Eq.(7) respectively, we obtain the equations 


of electric and magnetic field 


— BE, = (ViVi — V6) Ej, 
=g H= (VIV; — V6) Hy, (14) 


we read off the operator corresponding to L? from the above equations 


: VÍV; — V? 0 
D = Lf = , , (15) 
0 VÍV; — V? 
An eigenvalue of D is a nonnegative real number due to the Hermiticity of L. 
With Gauss law, we get two sets of equations taking the same form 
— OV; = DIV; = (VÍV: — V?) V, — V'Vi =, (16) 


where V; can be E; or H;. 


II. MIMICKING DE SITTER WITH METAMATERIALS 


Recently, the analogy between static curved space and inhomogeneous medium has been 
studied in A It is proved in [7| that in geometric optic limit, light propagates along the 
same path in electromagnetic material and in static curved spacetime as long as the optic 
metric 7;; is related to permittivity and permeability by Eq. (LJ). 

One can see the analogy more clearly as follows. To solve Maxwell equation in inhomo- 
geneous medium conveniently, the often used quantities are electric displacement field D! 


and magnetic induction field B’ defined as 
D’ = 6% Bi pe u” Hy, (17) 


where e” and u’ are permittivity and permeability specifying the electromagnetic properties 
of the medium. In terms of the newly defined quantities Dê and H*, Maxwell equations in 


inhomogeneous medium and Cartesian coordinate system are 
0:D'=0, ðB’ =0, (18) 
ðD — &8*0;H, =0, O,BY + &*0; Ey = 0. (19) 
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Since under spatial coordinate transformation €* and pi’ are tensor density with weight -1, 
we can define two tensors using €” and p by 
ij a ij 
E An o (20) 
For electromaganetic material with equal <% and p”, 77 = 7%. In this case, rewriting 
Eqs. (8) and in terms of 77, E; and H;, we find that Maxwell equations in the metame- 
terial takes the same form as Eq.) and Eq.(6) if 7’ is identified with the optic metric 
y 
By this analogy, one can design electromagnetic metamaterial corresponding a given met- 
ric. Along this line, theoretically, a new material called electromagnetic cloak was conceived 
36]. All parallel bundles of incident rays are bent around some region covered by electro- 
magnetic cloak and recombined in precisely the same direction as they entered the medium; 
Experimentally, a new class of electromagnetic materials are invented 8 which can be de- 
signed to have properties mimicking novel gravitational effects in Nature. Thus it is now 
conceivable that a material can be constructed whose permittivity and permeability values 
may be designed to vary arbitrarily. 
With the static de Sitter space 
2 ee 
ds? = —(1- Fae +(1- za) ar +r Po (21) 
The permittivity and permeability are 
06 90 sin 0 1 


rr rr Dai =Z = 
eT=p"=r'snd, E€” =p" =—,, il la T 
1- j (1 — =) sin 


(22) 
where (r,0, p) stands for spherical coordinate. In terms of Cartesian coordinate, the per- 
mittivity and permeability are 


5 eg 1 i¢ Tri 
ef =p = a0" - S. (23) 


The metamaterial can always be redesigned by making a coordinate mapping pi g. By 


making the following mapping 
r — f = Larcsin(r/L). (24) 


Then the permittivity and permeability becomes 


F 2Siml@/L) 00 sin 0 1 


8 00 — = PP — yP? — : 
cos(7/L) a j © j cos(ř/L) sin 8 


cos(f/L})’ (25) 


where (7,0, p) denotes the spherical coordinate. In terms of the Cartesian coordinate 


e = p = —__(§4 — = sin?(#/L) — 1) —_). (26) 


These two kinds of metamaterial mimicking the same de Sitter space, since they are 
related by a spatial coordinate transformation. However, the different permittivity and 
permeability will lead to distinct physical phenomenon in laboratory due to their different 
physical composition. 

The event horizon at r = L or fF = aL /2 now becomes the boundary of a cavity of 
metamaterial. As we shall see shortly, to metamaterial with € and pu as in Eq. (23), the leading 
term in the Casimir energy is a constant inversely proportional to certain microscopic length 
scale; while to metamaterial with € and u as in Eq.(26), the leading term in the Casimir 
energy in this cavity is proportional to L, we encourage experimenters to design such cavity 


to measure the Casimir energy, the result will have important implication for dark energy. 


IV. CASIMIR ENERGY OF ELECTROMAGNETIC FIELD IN STATIC DE SIT- 
TER SPACE 


In Maxwell theory, the electric sector and magnetic sector contribute the same amount 
to the Casimir energy, as the equations governing them are identical. In the following, we 
will focus on the electric sector. 


Consider the solution of electric field taking the form 
F(t, £) = e™E,(z), (27) 


where the absolute value of w should be interpretated as the energy with respect to time 
t. Upon quantization, the eigenvalue w contributes to the zero-point energy. The Casimir 


energy of electric field is the infinite sum 
E SD (28) 
asimir — F wW]. 
C 22 


Substitutting Eq.(27) into Eq. M6) we obtain the eigen equation of operator D with 
eigenvalue w? 


DIE; = wE; (29) 


Then in the framework of zeta function regularization, the problem of computing the Casimir 
energy is converted into the problem of rea heat kernel of D, the Casimir Energy 


can be extracted from the heat kernel as follows 


. 1 
E Casimir = lim z Ereg(+6) + Ereg(—€)}, (30) 
where 
1 
Hage 6(-5 +9) (31) 


with the zeta function associated with D defined by 


G(s) = a T ae ie). (32) 


In above expressions, a scale “jz” is introduced to keep the zeta function dimensionless; while 
the prime in tr’ indicates that zero eigenvalues are not included. In coordinate representa- 
tion, (a|e~*?|2’) is called the heat kernel denoted by K (t, x, x) satisfying the heat equation 
related to D as 

O,K(t,2,2)+DK(t,r,2) =0. (33) 


To solve this equation, we need to specify appropriate boundary conditions. 

At this moment, with above preparation, we can calculate the Casimir energy of electro- 
magnetic field in static de Sitter space time. This space time is described by the following 
metric 


2 2 
ds? = —(1— Fy) + — Fy) draen (34) 


where L is the de Sitter radius. However, the effective metric relevant to our calculation is 
the optic metric 4# appearing in the Maxwell equation with the following form 


r2 


2 
ds3 = yyde'de? = (1 — =) (1 - —) Pr + PO), (35) 
where the lower index “3” is used to emphasize that this metric is three dimensional. An 


interesting observation is that 7,; actually describes an anti de Sitter space. One can see 


this clearly by performing the following transformation on r coordinates: 


f 


r= —— (36) 


<2, 
ds? = (1 + dr + PEQ. (37) 


In terms of f, 


In the case of static de Sitter space time, the boundary condition Eq.(12) for retaining the 
Hermiticity of operator L in Eq.(i5) becomes 


(Fie Hog — Hio Eo¢)\r=r-5 = 9, (38) 


where the boundary is chosen at r = L— ô and 0 is a cut off to avoid the divergence of energy 
and its meaning will become clear later (see later for an explicit calculation of energy). Thus 


to keep Hermiticity, we choose 
Eo|r=L—s = Eplr=1— =0 (39) 


as the boundary condition. It is similar to the case of electromagnetic field in a spherical 
conductor shell. Therefore in the following, we will concentrate on the problem of finding 
the heat kernel corresponding to D defined with respect to yi; in Eq. (B5) and restricted to 


the space spanned by its eigenfunctions satisfying the Gauss law and boundary conditions. 


V'E; = 0, Pit = Eales = 0. (40) 


Usually, in a curved manifold the heat equation Eq. is hard to solve completely. However, 
the trace of the heat kernel has the following asymptotic expansion for small t, especially in 


three dimensions as 
i N 
tr (eP) = (—_)3/2 fow] bnt”) + olt”) }, 41 
P) = PPE eE fp t) (41) 


where a, are bulk terms independent of boundary condition, composed by polynomial of 
Riemann tensor, Ricci tensor, Ricci scalar, and their covariant derivatives; b,, are boundary 
terms comprised of the curvature, second fundamental form and their derivatives on the 
boundary. The sum runs over half-integers, but a, vanishes for half-odd-integers. Above 


asymptotic expansion implies that ¢(s) has a pole structure 


elo) = iy es + HS) (42) 


where f(s) is an entire analytic function of s, but in general, one has little information about 
it [9]. In the previous version of this work, we calculate the Casimir energy with above 
formulas. It is found that the leading term comes from a surface integral on the boundary 


horizon. Based on this expression, one may think that the Casimir energy is localized near 


the boundary horizon, but this may be a misunderstanding. For instance, in the frame work 
of heat kernel expansion, the Casimir energy between two infinite conductor plane in flat 
space will turn up as a boundary term depending on the permittivity and permeability of the 
conductor appearing to be localized on the boundary. However, L. Ford’s result T 
that in this case the Casimir energy distributes homogeneously between the two conductor 
planes. 

To avoid this misunderstanding, We will seak for a compact expression for the heat kernel. 
It is recalled that the optical metric describes an Euclidean AdS; space and the heat kernel 
defined with respect to it can be solved by utilizing Harmonic analysis |11]. To vector field 
the corresponding heat kernel satisfying the transverse condition and vanishing at the AdS3 
boundary was given by |12] (In our convention, the AdS; boundary corresponds to r = L 
or f = oo). From Eq. (40), it is obvious that the result of amounts to the 6 > 0 
limit of the heat kernel in our case. In other words, we can use the result of as a good 
approximation to the heat kernel sought by us up to the order 6/L. Based on that result 


we write down the contracted coincident limit of the heat kernel as 


2 


: t 2 
i K, (x, £, t) = 1+2—)e™/ 6/L). 43 
Y” Kylæ, x,t) an +273)€ " +o(ô/L) (43) 

With above result, tr (e~t?) can be evaluate by 
tr (e™tP) = | VTP Kalea ës 
L 2L 2 t 2 

= TL’(= — n(— 1+2—)e 1 44 
mL"(= a Gear Tage FoU) (44) 


where o(1) denotes terms < £. Substituting this result into Eqs. and (80), we obtain 


Fowinie = o—(in ye? = y — U(—5))(G — Fm) + of1/L) (45) 


where y is Euler constant and Y(—4) = I’ (—4)/T(-4). In above result, ô is a coordinate 
cut off so the physical meaning is not clear. To express it in terms of the physical cut off, 
we recall that in static de Sitter space the physical distance between r = L — ô and r = L 
is given by 

Alpny = L aresin(r/L)|"="_ 5. (46) 
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By requiring Alpny to be the Planck length lp, we derive 


l I" 
= L(1 — a eee 47 
§ = L(1 — cos(2)) = 5% (47) 
In terms of the physical cut off, the Casimir energy of electromagnetic field is 
3 1I L 1, 2L 
EGCasimir = — (l t= = = — = ln(— 1/L), 48 
caine = gD? = = VD gD) + 0/0 (48) 


where the dominant term proportional to L/ I? takes the same form as in the holographic 
dark energy model |2|. Compared with previous result given in lial kd, it is found that 
the leading term in the vacuum energy density of static de Sitter space is proportional to 
(L? —r?)~? which can also generate a leading term proportional to L/ iB after integrated over 
the bulk of static de Sitter space. Finally we note that the running renormalization In u? 


can be absorbed into k or be determined experimentally [13]. 


V. CASIMIR ENERGY OF ELECTROMAGNETIC FIELD IN ELECTROMAG- 
NETIC METAMATERIAL MIMICKING DE SITTER 


We have designed two kinds of metamaterial mimicking de Sitter in sect.3. For an elec- 
tromagnetic metamaterial with € and u taking the value of Eq. (@3) in laboratory, coordinate 
r has been the physical distance between origin and r, since metric measuring the physical 
distance is flat. In this case, the physical IR cut off is chosen to be L — d, and the Casimir 
energy is given by 


Ecim = z-a = 7 -YG -F a) + o(1/2), (49) 


where d is the counterpart of the Planck length in metamaterial, its value depends on the 
detailed chemical components and structure of the material which should be determined by 


experiment. We note that when L > d, 


ECasimir xy ir =y W(- )) i (50) 


1 
167 2 


al = 


L is usuaully much larger than d, the leading term is much larger than the usual term 
proportional to 1/L (for example, if d is nanometer, and L is 1cm, then the leading term is 


about 10° times of the second term, the usual term which is almost undetectable). 
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Alternatively, to metamaterial with € and pu taking the form as in Eq.(26), the Casimir 


energy is 


sin(ř/L) 1, ,1+sin(7/ZL) 


ECasimir T = mpri EL In(5 — sin(F/L) | 


(51) 
where 7 stands for physical distance between origin and 7 in this coordinate system and the 
size of the cavity is 7/2. So distinguished from the first case, 7L/2 — d is the physical IR 


cut off. d has the same meaning as in previous case. Usually L > d 
3 
E Casimir x zine ane io (=) ) 5- (52) 
T 


It is remarkable that in this metamaterial, the Casimir energy has the same form as in 


its gravity analogy, except that some microscopic scale d takes the place of Planck scale. 


VI. CONCLUSION 


It is a surprising result that the Casimir energy in a de Sitter space is not inversely 
proportional to the size of the horizon as for a usual finite cavity, rather it is proportional to 
the size of the horizon. However, if one is to expect that the Casimir energy is dark energy 
or at least a part of dark energy, this is a desired result, the normal answer in a cavity is 
just too small to be relevant in cosmology. 

Although we only computed the Casimir energy due to the photon field, it is expected 
that other fields will have the same form of Casimir energy. Since we suggest to construct 
a cavity of metamaterial to do experiment in laboratory, it is enough to know the form of 
the photon field Casimir energy, since this is what to be measured in laboratory. However, 
because the corresponding Casimir force is combined with the elastic force, it seems difficult 
to measure the predicted Casimir energy immediately. 

Aside from the issue of dark energy, to use metamaterials to mimic cosmology by itself 


is an exciting future direction, we look forward to developments in the future. 
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